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ABSTRACT 

Using Monte Carlo simulations, wc model the lum.inosity distribution of recycled pul- 
sars in globular clusters as the brighter, observable part of an intrinsic distribution 
and find that the observed luminosities can be reproduced using either log-normal or 
power-law distributions as the underlying luminosity function. For both distributions, 
a wide range of model parameters provide an acceptable match to the observed sample, 
with the log-normal function providing statistically better agreement in general than 
the power-law models. Moreover, the power-law models predict a parent population 
size that is a factor of between two and ten times higher than for the log-normal mod- 
els. We note that the log-normal luminosity distribution found for the normal pulsar 
population by Faucher-Giguere and Kaspi is consistent with the observed luminosi- 
ties of globular cluster pulsars. For Terzan 5, our simulations show that the sample 
of detectable radio pulsars, and the diffuse radio flux measurement, can be explained 
using the log-normal luminosity law with a parent population of ^ 150 pulsars. Mea- 
surements of diffuse gamma-ray fluxes for several clusters can be explained by both 
power-law and log-normal models, with the log-normal distributions again providing 
a better match in general. In contrast to previous studies, we do not find any strong 
evidence for a correlation between the number of pulsars inferred in globular clusters 
and globular cluster parameters including metallicity and stellar encounter rate. 

Key words: stars: neutron — pulsars: general — methods: numerical — methods: 
statistical — globular clusters: general — globular clusters: individual (Terzan 5) 
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1 INTRODUCTION 

The first millisecon d radio pulsar in a globular cluster (GC) 
was discovered by iLvne et al.l (jl987n . shortly after earlier 
predictions that the putative progenitors of millisecond pul - 
sars (MSP) are low-mass X-ray binaries (lAlpar et al.l 19821 ) 
which are known to be present in GCs ( Kata Il975l ) . In- 
spired by this discovery, a large number of sensitive pul- 
sar searches have been performed over the years resulting 
in the currently observed population of 143 radio pulsars 
in 27 GCiJ. GCs have some physical properties which are 
different than those of the galactic disk. Examples are ex- 
tremely high stellar density and high abundance of metal 
poor population II stars indicating that they were born in 
the early phase of the Galaxy's formation. These facts lead 
naturally to the question whether the population of radio 
pulsars in GCs is different to their counterparts in the Galac- 
tic disk. A number of differences are already well-known. 
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In particular, the high abundance of MSPs both in eccen- 
tric binary systems and as isolated objects. These phenom- 
ena can be explained as the results of two-body or three- 
body stellar interactions in the dense stellar environment of 
GCs jRomani. Kulk arni. fc Blandford"l987l: IVerbunt et al.l 
ll987l : llvanova et al.| l2008: Bagchi & Ray 20oi). On the other 
hand, the luminosity of a pulsar is a more fundamental prop- 
erty as it can in principle be linked to the pulsar emission 
mechanism. It is therefore important to establish whether 
there is any evidence for a different luminosity function of 
MSPs in the disk compared to those in GCs. 

There are two main ways to determine the pulsar lumi- 
nosity function numerically: (i) a full dynamical approach; 
(ii) a snapshot approach. In the dynamical approach, a sim- 
ulation is performed in which a model galaxy of pulsars is 
seeded according to various prescriptions of birth locations 
and initial rotational parameters. Each of these synthetic 
pulsars is then "evolved" both kinematically in a model for 
the Galactic gravitational potential and rotationally using 
a model for neutron star spin-down. The resulting popula- 
tion is then passed through the various detection criteria 
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(see, e.g.. lFaucher-Giguere fc KaspilbOOa '). The resulting set 
of "detectable" pulsars is then compared to the observed 
sample. The snapshot approach differs from the dynamical 
approach in that the pulsars are seeded at their final posi- 
tions in the galaxy without assuming anything about their 
spin-down or kinematic evolution and thus form a picture 
of the present day population. 

The dynamical approach has been used ex- 
tensive ly to study normal p ul sars in the g alactic 
disk |Bh attacharva et al. 199z; Gonthicr ct al.) [2002; 
iFaucher- G igucrc & Kaspi 2006; Ridlcv & Lorimorl [2010|). 
These studies considered the luminosity of pulsars to 
be described by power-law functions involving P and 
P with a substantial dispersion to account for distance 
un certainties and beam geometry. One of the conclusions 
by iFaucher-Giguere fc Kaspi ( 2006. hereafter FK06), also 
verified by iRidlev fc Lorimerl (|2010|) is that the resulting 
parent population of luminosities appears to be well de- 
scribed by a log-normal function. As discussed further in 
Section 3.1, the log-normal parameters favored by FK06 
(for the base- 10 logarithm of the 1400-MHz luminosities) 
are a mean of -1.1 and a standard deviation of 0.9. One of 
the goals of the current study is to examine whether these 
log-normal parameters are consistent with the observed 
populations of GC pulsars. 

For pulsars in GCs, where it is difficult to model the 
effects of stellar encounters and the cluster potential, the 
dynamical approach has so far not been used for radio pul- 
sar population syntheses. Although such an approach may 
be tractable in future, we will adopt in this work a version of 
the snapshot approach. We will carry out Monte Carlo sim- 
ulations that assume all GCs have the same intrinsic pulsar 
luminosity function, but a different population size and use 
this approach to model the observed sample of pulsars given 
the various ranges of luminosities. As we shall see, this ap- 
proach provides a remarkably good agreement between the 
model and observed luminosity distributions. 

Within the snapshot framework, one usually fits the 
complementary cumulative distribution function (CCDF) of 
pulsar luminosity as a power law as A'^(> I/1400) = NoL\^qq, 
where Z/1400 is the luminosity of the pulsar at 1400 MHz, 
N is the number of pulsars having luminosity value greater 
than 1/1400, Nq and q are constants. We note here that some- 
times in literature (JHessels et al.ll2007l : iHui. Cheng fc TaamI 
[2OI0I), the CCDF has been mentioned as the cumulative 
distribution function; but actually the cumulative distri- 
bution function (CDF, N{^ ii4oo) ) is related to the 
CCDF as CDF = 1- CCDF . Using this snapshot technique, 
iHui. Cheng fc TaamI l|2010f ) concluded that the luminosities 
of MSPs in GCs are different from those in the galactic disk 
as the CCDF for GC MSPs is much steeper than that of 
disk pulsars. This is a very important conclusion. If correct, 
it would imply that the radio luminosity is related to dif- 
ferences in formation processes between the disk and GC 
pulsars. The same analysis was re-performed with more re- 
cent distance estimates of GCs and the resultant CCDF was 
even steeper (Bagchi fc Lorimer 2010). 

Since GCs are generally at large distances, the lumi- 
nosity function of observed pulsars is not as well sampled 
as in the Galactic disk. In the present work we try to ac- 
count for this incompleteness by considering GC MSPs as 
the brighter tail of some intrinsic parent population. The 



goal of the current work is to explore the range of possible 
distributions that are consistent with the current sample of 
GC pulsars. The plan for the rest of this paper is as fol- 
lows. In Section [2I we describe the pulsar sample we use. In 
Section [31 we present our analysis procedure and its main 
results. In Section[3| we investigate additional constraints on 
allowed model parameters from observations of diffuse radio 
and gamma-ray flux. In Section [SI we compare our results 
with earlier work. We draw our main conclusions in Section 

m 



2 THE OBSERVATIONAL SAMPLE 

The luminosity of a puls ar, L, can be computed (see, e.g., 
[Lorimer fc Krameij [2005[ ) from its distance and the mean 
flux density S^ (defined at some observing frequency v) us- 
ing the following geometrical relationship: 



And'' 



Svdv. 



(1) 



Here p is the radius of the (assumed circular) emission cone, 
5 is the pulse duty cycle (= Weq/P), P is the spin period 
of the pulsar and Weq is the equivalent width of the pulse 
(i.e. the width of a top-hat shaped pulse having the same 
area and peak flux density as the true proflle), uiovi and v^p 
deflnes the range of radio frequencies over which the pulsar 
has been observed and d is the distance of the pulsar. As it is 
usually difficult to determine the values of p and 5 reliably, 
we define the "pseudoluminosity" 



Li^ — o^ d 



(2) 



Henceforth, we will use the term luminosity to mean pseu- 
doluminosity. 

Among the current sample of 143 GC pulsars, flux den- 
sity values have been reported for 107 pulsars. Among these, 
three are clearly young isolated objects which more closely 
resemble the normal population of pulsars i n the Galaxy 
(for a further discussion of this population, see [Bovles et al.[ 
[2OIH ). We consider here the sample of 83 pulsars in 10 GCs 
with spin periods P ^ 100 ms and each of these GCs host at 
least 4 such pulsars. For all these objects, the spin and bi- 
nary properties suggest that the neutron star has undergone 
a phase of recycling in the past. 

Among these, for 45 pulsars (14 in 47 Tuc, 4 in M3, 
5 in M 5, 5 in M 13, 5 in NGC 6752, 3 in NGC 6517 and 
9 in M 28) flux density values have been measured at 1400 
MHz; for 31 pulsars (25 in Terzan 5, 5 in NGC 6440 and 1 in 
NGC 6517) flux density values have been measured at 1950 
MHz and for 7 pulsars (in M 5) flux density values have been 
measured at 400 MHz. To pursue our study of pulsar lumi- 
nosities at 1400 MHz, we scale the flux densities measured 
at other frequencies using the power-law Si, oc v" , where 
Q is the spectral index. We then use the model prediction 
for 51400 as the best estimate of the pulsar's flux density at 
this frequency. In these calculations we use the estimated 
values of a from observed values of fluxes at different fre- 
quencies whenever available, otherwise adopt the mean a 
of GC MSPs (for which flux values have been rep orted at 
multiple frequencies) of — 1.9. [Toscano et al.[ ([l998[ ') also ob- 
tained mean a of 19 millisecond pulsars to be —1.9 ± 0.1, 
but their sample contains only two GC pulsars. Once we get 
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Table 1. Fluxes and spectral indices of 107 pulsars in globular clusters, a can be calculated for 20 pulsars using the central frequency of 
observations (when ever reported). Mean a of pulsars having Ps ^ 100 ms is -1.865 (excluding positive a of PSR J1836— 2354A). We set 
a = —1.9 in the present work. From the sample of 107 pulsars, we exclude pulsars with spin period > 100 ms; and then exclude pulsars 
for which the host GC contains less than 4 pulsars with Ps ^ 100 ms and known flux values. Pulsars which are not used in the present 
study have been written in italics. References are at the end of the table. 



GC 


PSR 


Ps 


5400 


•Seoo 


51170 


51400 


•Sieoo 


5'l950 


a 






(ms) 


(mjy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 




47Tuc 


J0023-7204C 


5.757 


1.53(ri) 


1.54(rJ) 


* 


0.36{r2) 


* 


* 


-1.352 


47Tuc 


J0024-7204D 


5.358 


0.95(ri) 


0.55(r;) 


* 


0.22{r2) 


* 


* 


-1.264 


47 Tuc 


J0024-7205E 


3.536 


* 


* 


* 


0.21(r2) 


* 


* 


* 


47Tuc 


J0024-7204F 


2.624 


* 


* 


* 


0.15(r2) 


* 


* 


* 


47 Tuc 


J0024-7204G 


4.040 


* 


* 


* 


0.05(r2) 


* 


* 


* 


47 Tuc 


J0024-7204H 


3.210 


* 


* 


* 


0.09(r2) 


* 


* 


* 


47 Tuc 


J0024- 72041 


3.485 


* 


* 


* 


0.09(r2) 


* 


* 


* 


47 Tuc 


J0023-7203J 


2.101 


* 


* 


* 


0.54(r2) 


* 


* 


* 


47 Tuc 


J0024-7204L 


4.346 


* 


* 


* 


0.04(r2) 


* 


* 


* 


47 Tuc 


J0023-7205M 


3.677 


* 


* 


* 


0.07(r2) 


* 


* 


* 


47 Tuc 


J0024-7204N 


3.054 


* 


* 


* 


0.03(r2) 


* 


* 


* 


47 Tuc 


J0024-7204O 


2.643 


* 


* 


* 


0.10(r2) 


* 


* 


* 


47 Tuc 


J0024-7204Q 


4.033 


* 


* 


* 


0.05(r2) 


* 


* 


* 


47 Tuc 


J0024-7203U 


4.343 


* 


* 


* 


0.06(r2) 


* 


* 


* 


NGC 1851 


J0514-4002 


4.990 


0.28(r5) 


* 


* 


* 


* 


0.0056(r5) 


-2.568 


M 53 


B1310+18 


33.163 


1.0{r4) 


* 


* 


* 


* 


* 


* 


M3 


J1342+2822A 


2.545 


* 


* 


* 


0.007(r5) 


* 


* 


* 


M3 


J1342+2822B 


2.389 


* 


* 


* 


0.014(r5) 


* 


* 


* 


M3 


J1342+2822C 


2.166 


* 


* 


* 


0.006(r5) 


* 


* 


* 


M3 


J1342+2822D 


5.443 


* 


* 


* 


0.010(r5) 


* 


* 


* 


M 5 


B1516+02A 


5.554 


0.5(r6) 


* 


0.155(r7) 


0.120(r5) 


* 


* 


-1.161 


M 5 


B1516+02B 


7.947 


0.3{r6) 


* 


0.027(r7) 


0.025(r5) 


* 


* 


-2.132 


M 5 


J1518+0204C 


2.484 


* 


* 


* 


0.039(r5) 


* 


* 


* 


M 5 


J1518+0204D 


2.988 


* 


* 


* 


0.008(r5) 


* 


* 


* 


M 5 


J1518+0204E 


3.182 


* 


* 


* 


0.010(r5) 


* 


* 


* 


M 4 


B 1620-26 


11.076 


15(rS) 


7.2(r5) 


* 


1.6{rlO) 


* 


* 


-1.744 


M 13 


B1639+36A 


10.378 


3.0(r^) 


* 


* 


0.140(r5) 


* 


* 


-2.486 


M 13 


B1639+36B 


3.528 


* 


* 


* 


0.022(r5) 


* 


* 


* 


M 13 


J1641+3627C 


3.722 


* 


* 


* 


0.030(r5) 


* 


* 


* 


M 13 


J1641+3627D 


3.118 


* 


* 


* 


0.024(r5) 


* 


* 


* 


M 13 


J1641+3627E 


2.487 


* 


* 


* 


0.010(r5) 


* 


* 


* 


M 62 


J1701-3006A 


5.242 


* 


* 


* 


OA{rll) 


* 


* 


* 


M 62 


J1701-3006B 


3.594 


* 


* 


* 


0.3(rll) 


* 


* 


* 


M 62 


J1701-3006C 


7.613 


* 


* 


* 


0.3{rll) 


* 


* 


* 


NGC 6342 


B1718-19 


1004.04 


0.2b3{rl2) 


0.550(riS) 


* 


0.278(ri2) 


0.1S{rl2) 


* 


-0.338 


NGC 6397 


J 1740-5340 


3.650 


* 


* 


* 


l.0{rl3) 


* 


* 


* 


Ter 5 


J1748-2446A 


11.563 


* 


5{rl4) 


* 


0.61{rl5) 


* 


1.020{rl6) 


-1.572 


Ter 5 


J1748-2446C 


8.436 


* 


* 


* 


* 


* 


0.360{rl6) 


* 


Ter 5 


J1748-2446D 


4.714 


* 


* 


* 


* 


* 


0.041 {rl6) 


* 


Ter 5 


J1748-2446E 


2.198 


* 


* 


* 


* 


* 


0MS{rl6) 


* 


Ter 5 


J1748-2446F 


5.540 


* 


* 


* 


* 


* 


0m5{rl6) 


* 


Ter 5 


J1748-2446G 


21.672 


* 


* 


* 


* 


* 


0.015{rl6) 


* 


Ter 5 


J1748-2446H 


4.926 


* 


* 


* 


* 


* 


0.016{rl6) 


* 


Ter 5 


J1748-2446I 


9.570 


* 


* 


* 


* 


* 


0.029{rl6) 


* 


Ter 5 


J1748-2446J 


80.338 


* 


* 


* 


* 


* 


0.019{rl6) 


* 


Ter 5 


J1748-2446K 


2.970 


* 


* 


* 


* 


* 


0.040{rl6) 


* 


Ter 5 


J1748-2446L 


2.245 


* 


* 


* 


* 


* 


0.041{rl6) 


* 


Ter 5 


J1748-2446M 


3.570 


* 


* 


* 


* 


* 


0m3{rl6) 


* 


Ter 5 


J1748-2446N 


8.667 


* 


* 


* 


* 


* 


0m6{rl6) 


* 


Ter 5 


J1748-24460 


1.677 


* 


* 


* 


* 


* 


0.120{rl6) 


* 
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Table 1. (continued). 



GC 


PSR 


Ps 


S'400 


'S'eoo 


5'ii7a 


5'i400 


'S'leoo 


Sigso 


a 






(ms) 


(mjy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 




Ter 5 


J1748-2446P 


1.729 


* 


* 


* 


* 


* 


0.077{rJ6) 


* 


Ter 5 


J1748-2446Q 


2.812 


* 


* 


* 


* 


* 


0.027{rJ6) 


* 


Ter 5 


J1748-2446R 


5.028 


* 


* 


* 


* 


* 


0.012{rJ6) 


* 


Ter 5 


J1748-2446S 


6.117 


* 


* 


* 


* 


* 


0.018{rJ6) 


* 


Ter 5 


J1748-2446T 


7.085 


* 


* 


* 


* 


* 


0.020{rJ6) 


* 


Ter 5 


J1748-2446U 


3.289 


* 


* 


* 


* 


* 


0.016{rJ6) 


* 


Ter 5 


J1748-2446V 


2.072 


* 


* 


* 


* 


* 


0.071(ri6) 


* 


Ter 5 


J1748-2446W 


4.205 


* 


* 


* 


* 


* 


0.022(ri6) 


* 


Ter 5 


J1748-2446X 


2.999 


* 


* 


* 


* 


* 


0.018(ri6) 


* 


Ter 5 


J1748-2446Y 


2.048 


* 


* 


* 


* 


* 


0.016{rJ6) 


* 


Ter 5 


J1748-2446ad 


1.396 


* 


* 


* 


* 


* 


0.08{rJ7) 


* 


NGC 6440 


B 1745-20 


288.603 


10{rl8) 


* 


* 


0.37{rl5) 


1.5{rJS) 


0.37{rJ9) 


-1.920 


NGC 6440 


J1748-2021B 


16.760 


* 


* 


* 


* 


* 


0.047{rJ9) 


* 


NGC 6440 


J1748-2021C 


6.227 


* 


* 


* 


* 


* 


0.044{rJ9) 


* 


NGC 6440 


J1748-2021D 


13.496 


* 


* 


* 


* 


* 


0.075{rJ9) 


* 


NGC 6440 


J1748-2021E 


16.264 


* 


* 


* 


* 


* 


0.023{rJ9) 


* 


NGC 6440 


J1748-2021F 


3.794 


* 


* 


* 


* 


* 


0.017{rJ9) 


* 


NGC 6441 


J 1750- 37 A 


111.608 


* 


* 


* 


* 


* 


0.059(rJ9) 


* 


NGC 6441 


J 17 50- 3703 B 


6.074 


* 


* 


* 


* 


* 


0.037(rJ9) 


* 


NGC 6441 


J1750-3703C 


26.569 


* 


* 


* 


* 


* 


0.015(rJ9) 


* 


NGC 6441 


J 17 50- 3703 D 


5.140 


* 


* 


* 


* 


* 


0.010{rJ9) 


* 


NGC 6517 


J1801-0857A 


7.176 


* 


* 


* 


0.036(r20) 


* 


0.020{r.80) 


-1.648 


NGC 6517 


J1801-0857B 


28.961 


* 


* 


* 


0.012(r20) 


* 


0.009(raO) 


-0.806 


NGC 6517 


J1801-0857C 


3.739 


* 


* 


* 


0.012(r20) 


* 


0.007(r.80) 


-1.511 


NGC 6517 


J1801-0857D 


4.226 


* 


* 


* 


* 


* 


0.011(r.gO) 


* 


NGC 6539 


B1802-07 


23.101 


3.1{r21) 


1.0{r21) 


* 


0.6 {r21 ) 


* 


* 


-1.213 


NGC 6544 


J1807-2459A 


3.059 


* 


* 


* 


l.3ir22) 


* 


* 


* 


NGC 6624 


B1820-30A 


5.440 


W{r9) 


6.8{r9) 


* 


0.72{r9) 


0.31(r25) 


* 


-2.922 


NGC 6624 


B1820-30B 


378.596 


2.2{r21) 


1.0{r21) 


* 


0.07(r25) 


0.07(r25) 


* 


-2.654 


M 28 


B1821-24A 


3.054 


30ir24) 


* 


* 


0.94(r25) 


* 


* 


-2.764 


M 28 


J1824-2452B 


6.547 


* 


* 


* 


0.07(r25) 


* 


* 


* 


M 28 


J1824-2452C 


4.159 


* 


* 


* 


0.17(r25) 


* 


* 


* 


M 28 


J1824-2452D 


79.832 


* 


* 


* 


0.05(r25) 


* 


* 


* 


M 28 


J1824-2452E 


5.420 


* 


* 


* 


0.06(r25) 


* 


* 


* 


M 28 


J1824-2452F 


2.451 


* 


* 


* 


0.08(r25) 


* 


* 


* 


M 28 


J1824-2452G 


5.909 


* 


* 


* 


0.05(r25) 


* 


* 


* 


M 28 


J1824-2452H 


4.629 


* 


* 


* 


0.06(r25) 


* 


* 


* 


M 28 


J1824-2452J 


4.039 


* 


* 


* 


0.07(r25) 


* 


* 


* 


M 22 


J1836-2354A 


3.354 


* 


* 


* 


0.040(r20) 


* 


0.043(ra0) 


0.203 


M 22 


J1836-2354B 


3.232 


* 


* 


* 


0.200(r20) 


* 


0.073(r.gO) 


-2.826 


NGC 6749 


J1905+0154A 


3.193 


* 


* 


* 


0.023(r5) 


* 


* 


* 


NGC 6749 


J1905+0154B 


4.968 


* 


* 


* 


0.006(r5) 


* 


* 


* 


NGC 6752 


J1911-5958A 


3.266 


* 


* 


* 


0.21(r26) 


* 


* 


* 


NGC 6752 


J1910-5959B 


8.358 


* 


* 


* 


0.05(r26) 


* 


* 


* 


NGC 6752 


J1911-6000C 


5.277 


* 


* 


* 


0.24(r26) 


* 


* 


* 


NGC 6752 


J1910-5959D 


9.035 


* 


* 


* 


0.05(r26) 


* 


* 


* 


NGC 6752 


J1910-5959E 


4.572 


* 


* 


* 


0.07(r26) 


* 


* 


* 


M 71 


J 1953 -1-1846 A 


4.888 


* 


* 


* 


0.059(r5) 


* 


* 


* 


M 15 


B2127-hllA 


110.665 


1.7{r27) 


* 


* 


0.2{r.8S) 


* 


* 


-1.797 


M 15 


B2127+11B 


56.133 


1.0{r27) 


* 


* 


* 


* 


* 


* 
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Table 1. (continued). 



GC 



PSR 



Ps 

(ms) 



•S'400 
(mjy) 



■^600 ■S'1170 

(mJy) (mJy) 



•S'l400 
(mJy) 



■51600 •S'lgso 
(mJy) (mJy) 



M 15 


B2127+11C 


30.529 


0.64(r27) 


M 15 


B2127+11D 


4.803 


0.34(r27) 


M 15 


B2127+11E 


4.651 


0.24(r27) 


M 15 


B2127+11F 


4.027 


0.14(r27) 


M 15 


B2127+11G 


37.660 


0.13(r27) 


M 15 


B2127+11H 


6.743 


0.16(r27) 


M 30 


J2140-2310A 


11.019 


* 



0.08(r29) 



References : (1) rl i lRobinson et al.l lll995l). ( 2) r2 i lCamilo et"all ||2000|'), (3) rS : iFreire. Ransom 



1)^ 



iKulkarni et al.l lll99li '). (5)^5 : lHessels et al.l ll2007tl . (6) r6 : lAnderson et al.1 19971). (J) r7 ilFreire et al. ll2008al). (8) r8 : un published 
( http://www. atnf. csiro . au/research/pulsar/psrcat/expert.html), (9) r9 : lToscano et al.l 1119981). (10) r lO : Kramer et al.l 11 19981). ( 11 



Guptal 1I2OO7I) . (4) r4 



: |Possenti et al.l ll2003l), (12) rl2 : Aver aged over v ariations with the o rbital phase llLvne et al.lil993l ). ( 13) rl3 :lD'Amico et al.1 |20o3) 
(14) rU :lLvne et alj lll990l ), (15) rl5 : Hobbs ct al. (20gj), (16) rl6 ilRansom et aL ||2005'). (17) rl7 ilHessels et al.' ('2006"). (18) rl8 

iLvne et al.llll996l) . (19) rl9 JPrcirc ct al. (2008b), (20) r2 :lLvnch et al.l 1I2OIM. (21) r21 : lLorimer et al.l lll 995) 
iRansom et al.l 1I2OOII). (23) r 23 i lsiggs et al.l lll994|) (24) r24 : lFoster. Fairhead. fc Backed lll99ll ) 
ICorongiu etldl tOOdh . (27) r27 i lAndersonI lll993l) . (28) r28 : IWolszczan et al.l lll989l) 



(22) r22 : 

(25) r25 :lBeginlli2 006). (26) r26 
(29) r29 : iRansom et al., (2004 ) 



100 



10 : 





a=-1.9 


^ "v.,^ 








; + 






^ 



Li 400 ['"Jy kpc 1 

Figure 1. Complementary cumulative distribution of 1400-MHz 
luminosities for the sample of globular cluster pulsars. Flux den- 
sity values measured at other frequencies have been converted 
to S1400 using Si, oc u" for a = —1.9. We have checked that a 
change in the value of « ~ ±0.3 does not make any visible change 
in the shape of the plot. 



5i4oo, we can calculate I/1400 if d is known. For GC pulsars, 
the distances are taken to be those of their host clusters. We 
use the most recent distance estimates from the literature. 
Our complete list of GC pulsar flux and spectral param- 
eters used in this section, and the remainder of the paper is 
given in Table[T] While compiling this list, we confirmed the 
earlier conclusions by Hesscls ct al. (2007) that the choice 
of a in a realistic range does not affect the complementary 
cumulative distribution (CCD) of luminositie s signif icantly. 
Using 37 isolated GC pulsars, iHessels et al.l (|2007l ) found 
that an arbitrary choice of a in the range of -1.6 to -2.0 does 
not affect the shape of the CCD. We arrive at the same con- 
clusion with our sample using different values of a as -1.6, 
-1.9 and -2.2. We perform Kolmogorov-Smirnov (KS) tests 
between the luminosity distributions obtained with different 
choice of a. The KS test returns a statistic Pks which gives 
the probability that the two samples are drawn from the 



same distribution (for details, see lPress et al. 1120071 ). In this 
case, Pks is always greater than 0.997 when we compare any 
two luminosity distributions from the three obtained with 
Q = —1.6, —1.9 and —2.0; as an example, the distribution 
obtained with a = —1.9 is shown in Fig. [l] 

The absence of any significant difference in the lu- 
minosity distribution for a realistic range of a also 
supports our choice of a = —1.9 for this work, 
which is not too different than the favor ed choice 
(a = -1.8) of other studies (iHui. Cheng fcTa am 20l3; 
iMaron. Kiiak. Kramer fc Wielebinskil 119871 ). As we are 



studying luminosities only at 1400 MHz, henceforth, we shall 
denote Z/1400 simply by L. 

Note also that, throughout this paper, we will be con- 
cerned with the possible forms of the luminosity distribution 
of cluster pulsars. Correlations between luminosity and other 
pulsar parameters are not discussed in any detail here. The 
reason for this is that, as for the pulsar population in the 
Galaxy correlations in the observed pulsar samples are not 
apparent due to the presence of distance errors and beam- 
ing uncertainties (see Lorimer et al. 1993 for a discussion). 
We did not find any correlation between luminosity and spin 
period of the 83 recycled pulsars used in the present work. 
As mentioned earlier, P measurements for globular cluster 
pulsars are affected by cluster potential, so can not easily 
be used to study intrinsic properties of the pulsars. For the 
remainder of this paper, we proceed with the underlying as- 
sumption that there exists a single luminosity function for 
all globular cluster pulsars, and attempt to explain the ob- 
served luminosities in this way. As we will demonstrate, the 
data are remarkably consistent with this simple idea. How- 
ever, the wide ranges of possible model parameters that are 
consistent with the data do not rule out the idea that the 
parent luminosity function may vary from cluster to cluster. 



3 ANALYSIS 

With the data described above, we aim to find luminosity 
distribution functions whose brighter tail can be considered 
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as the observed luminosity distribution of GC pulsars, as- 
suming that the parent luminosity distribution is the same 
for all GCs. To do so, for each GC, we first generate a syn- 
thetic sample of A^'triai.i pulsar luminosities from a chosen 
distribution function until we get C x Nohs,i pulsars having 
simulated luminosities greater than the observed minimum 
luminosity for that GC. This multiplication by the constant 
'C (100-1000) is done to minimize statistical variations. In 
this notation, i is the GC index NahB,i is the observed num- 
ber of pulsars in the GC that we consider. 



C 



E 



(3) 




is the total luminosity and 



-'sim.tot — 



C 



/ , '-'simj 



is the total flux in the i' GC. Here Ls_ 



(4) 



and Ssim,i are 



the simulated luminosities and corresponding fluxes. After 
we perform the simulation for all 10 GCs, we compare the 
simulated luminosities with the observed luminosities of 83 
pusars by performing KS and y^ tests. As mentioned earlier, 
the KS test can be used to test the hypothesis that two 
distributions differ, with a low value of KS probability Pks 
suggesting a mismatch. The x^ statistic uses binned data 
and compares the values of the two distributions at each 
bin; here a low value of y^ implies a good agreement. Here 
we divide the luminosity range 0.1-1000 mjy kpc^ into 36 
logarithmically equispaced bins. A'triai,i/C is the predicted 
number of total pulsars in that GC which we call as A^rad,!. 
A key assumption in our present analysis is that each 
GC has been searched down to the level of the faintest ob- 
servable pulsar in that particular cluster. This assumption 
provides a good approximation to the actual survey sensi- 
tivity in each cluster, and was made primarily due to the 
lack of currently published detail of several of the globular 
cluster surveys so far. The assumption greatly simplifies our 
modeling procedure, since it means that we do not have to 
consider variations in sensitivity due to other factors (for 
example scintillation, eclipsing binary systems etc.). This 
simple approach is appropriate for the purposes of the cur- 
rent work where we are simply trying to assess the range 
of luminosity functions compatible with the data. A more 
rigorous study which takes account of the survey thresholds 
in detail may well be able to narrow the range of possible 
model parameters found here, and should certainly be car- 
ried out when more details of the surveys are published, but 
is beyond the scope of the current work. 

3.1 Log-normal luminosity function 

We begin by testing a log-normal luminosity function, where 
the probability density function (PDF) 



h 



log — norm. 



alW 



log 



^ 



exp 



-(log 



2a2 



(5) 

where, as usual, fj, is the mean of the distribution and a is 
the standard deviation. For this choice of distribution, we 
find that C = 100 is sufficient to minimize statistical fluctu- 
ations. The variation of Pks and x^ with ^ and a are shown 




Figure 2. Top: variation of Pks (for Pks ^ 0.05) with ^ and a. 
Bottom; variation of x^ (within 2u about the minimum value of 
X^) with /i and a, the parameters of the log-normal distribution. 
A 1(T contour is also shown. 



in Fig. (2] It is clear that there is a wide range of values of 
yL, a for which the simulated luminosity distributions agree 
well with the observed sample. For the two statistical tests, 
good agreement is given when Pks has high values and x'^ 
is small. As expected, the region oi [i — a parameter space 
encompassed by Pks > 0.05 is essentially the same as the 
contours encompassing the 95% probability values around 
the x^ minimum. 

For this distribution, and for the purposes of later dis- 
cussion, we deflne three models based on particular param- 
eter choices. Model 1 uses the parameters found by FK06 
(/i = —1.1 and a = 0.9) from which we flnd Pks = 0.15 
and x^ = 9.4. Model 2, for which p. = -0.61 and a = 0.65 
returns the maximum value of Pks = 0.98 with a x^ ~ 7.9. 
Model 3, for which /i — —0.52 and a — 0.62, returns a min- 
imum value of x^ ~ 6.3 and has Pks = 0.37. 

In Fig. [3] we compare these three models with the ob- 
served data. As expected, all models match well. While 
model 3 provides the closest match by eye, the statistical 
results mentioned above do not rule out either model 1 or 
model 2. The FK06 luminosity model parameters (model 1), 
therefore, are consistent with the observed CCD. 



3.2 Power-law luminosity function 

As mentioned earlier, power-law luminosity functions have 
been used by a number of authors. It is therefore of great 
interest to see how the power-law compares to log-normal 
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L [mjy kpc"] 




L [mJy kpc"] 



Figure 3. Observed luminosity distribution witli simulated lumi- 
nosity distributions generated with log-normal distribution func- 
tion for three different sets of fi and a which are defined as models 
1 (upper curve), 2 (lower curve) and 3 (middle curve); see text 
for details. "*"s represent the observed distribution. 



for the GC pulsars. The PDF of the power law distribution 



/po 






(6) 



where Lmin is the minimum value of L and /3 the power-law 
index. This abrupt cut-off required to avoid divergence when 
integrating this function over all L is somewhat unphysical, 
but nevertheless can be used to parameterize the luminosi- 
ties in an independent way to the log-normal. We perform 
simulations over a range of I/min 0.003 — 0.48 mJy kpc , as 
0.48 mJy kpc^ is the observed minimum luminosity among 
GC pulsars in our sample, and the lower value of Lmin is 
chosen somewhat arbitrarily. For this model, we found that 
C = 1000 is required to minimize statistical fluctuations. 
Unlike the log-normal model, we found that the power law 
distribution occasionally produced pulsars with large lumi- 
nosities L ^ 100 mJy kpc'^ which biased some of our prelim- 
inary simulation runs. To avoid this difficulty, we imposed 
a maximum luminosity of 50 mJy kpc^. No GC pulsar is 
currently known with L > 20 mJy kpc^, and our results are 
insensitive to the exact choice of the maximum luminosity 
cutoff over the range 20-500 mJy kpc^. 

The nominal best parameter values give Pks = 0.81 
{X^ = 8.7) for l3 = 0.92, Lmin = 0.017 mJy kpc^ (model 
4) and minimum x^ = 8.0 (Pks = 0.56) for /3 = 1.01, 
imin = 0.022 mJy kpc^ (model 5). Our values of /3 (which 
give good fits) are not too different from the conventional 
values P -| - 1 = 2. For example, the best fit value of the 
analysis of iFruchter fc Gosa (2000) for Terzan 5 pulsars was 
fi + 1 = 1.85. We have seen that for /3 ~ 1, the fit does 
not dep end much on Lm i n whic h again agrees with the re- 
sults of IFruchter fc Gosa l|2000f ). As an additional point of 
reference, we als o consider the nominal p ower-law parame- 
ters discussed by IFruchter fc GossI (|2000| ). i.e. /3 = 0.85 and 
Lmin = 0.03 mJy kpc^ which we refer to as model 6. For 
this pair ing, Pks = 0.43 and y^ = 9.0. We need to remem- 
ber here IFruchter fc Gosd (|2000|) did not put any constraint 
on the maximum value of the luminosity. In Fig. |3] we com- 
pare the observed CCD with simulated CCDs for models 
4-6. Statistically, the agreement between simulated and ob- 



Figure 4. Observed luminosity distribution with simulated lu- 
minosity distributions generated with a power-law distribution 
function for three different set of Lmin and /3 which are defined 
as models 4 (middle curve), 5 (lower curve) and 6 (upper curve); 
see text for details. "*"s represent the observed distribution. The 
mismatch between the simulated and observed CCDs appear only 
when the number of pulsars are < 10. 



served distributions is almost as good as that for log-normal 
distributions. 



3.3 Exponential Distribution 

The above models characterized the luminosity function in 
terms of two parameters. For completeness, we also consider 
a simple one-parameter model, the exponential distribution 
with PDF 



/. 



exponential ^-^^1400 



(Li 



Ae" 



(7) 



Here 1/A is the mean of the distribution, and we find that 
C = 100 is enough to get rid of statistical fluctuations. 
For this model, the maximum value of Pks = 0.17 is ob- 
tained for A = 0.676 mjy~^ kpc~^ with a corresponding 
X^ ~ 13.7. The minimum value of x^ ~ 11-94 is found for 
A = 0.439 mJy"^ kpc"^ with a Pks = 0.00028. In Fig.[5]we 
compare the observed CCD with the simulated CCDs with 
these two values of A. It is clear that the simulated distri- 
bution never matches with the observed one very well. We 
therefore do not consider the exponential distribution fur- 
ther in this work and focus the remainder of the discussion 
on the log-normal and power-law distributions. 



4 MODEL PREDICTIONS AND 
CONSTRAINTS 

In the previous section, we found that there is a large family 
of possible luminosity parameters that are consistent with 
the observed distribution of GC pulsar luminosities. These 
ranges translate to a variety of different predictions for the 
population sizes in each GC. This can be seen from a com- 
parison of the predicted parameters for each GC using the 
log- normal parameter choices (models 1, 2 and 3) in Table [2] 
and the power-law parameter combinations (models 4, 5 and 
6) in Table O In this section, we try to place further con- 
straints on these parameters by examining the predictions 
for the diffuse radio and gamma-ray fluxes separately. 
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Table 2. Population estimations and predictions using log-normal luminosity functions. For each cluster, we list the predicted number of 
potentially observable radio pulsars {Af^ad): the predicted total diffuse radio flux (Saim.tot)! the predicted gamma-ray luminosity (i^,sim) 
for three different choices of efficiency (177). See text for further details. To compute uncertainties in A^radi wo assume that they are 
dominated by the statistical noise in the observed number of pulsars, A^obs- The uncertainty in N^^^ is then simply N^^^/^JNahs- Also 
listed for each model is A'^10, the total population estimate for these 10 GCs. 



Cluster 


Af^ad 


'-'sim,tot 

(mjy) 


{■q^) = 0.08 


-^7,sim 

(r,7> = 0.06 
(1034 erg s-i) 




(r,7>=0.1 




Model 1 (FK06): ^ 


= — 1.1 and a = 


0.9 A^io = 688 ± 82 




47Tuc 


71 ±19 


3.1 ±0.8 


10 ±5 


7.7 ±3.6 




13 ±6 


M3 


24 ± 12 


0.16 ±0.08 


3.5 ±2.2 


2.6 ±1.6 




4.3 ±2.7 


M5 


24 ±11 


0.31 ±0.14 


3.5 ±2.1 


2.6 ±1.6 




4.3 ±2.6 


M13 


25 ±11 


0.38 ±0.16 


3.6 ±2.1 


2.7 ±1.6 




4.5 ±2.6 


Ter 5 


167 ±33 


3.7 ±0.7 


24 ± 11 


18 ±8 




30 ±13 


NGC 6440 


88 ±39 


0.86 ±0.4 


13 ±8 


10 ±6 




16 ±9 


NGC 6517 


46 ±23 


0.29 ±0.15 


6.6 ±4.2 


5.0 ±3.1 




8.3 ±5.2 


M28 


120 ±40 


2.6 ±0.9 


17 ±9 


13 ±7 




22 ±11 


NGC 6752 


44 ±20 


1.7 ±0.8 


6.3 ±3.8 


4.8 ±2.8 




7.9 ±4.7 


M15 


79 ±30 


0.52 ±0.20 


11 ±6 


8.5 ±4.6 




14 ±8 




Model 2 ( 


[maximum P^s) 


: 11 = —0.61 and 


a = 0.65 Afio = 453 ± 56 




47Tuc 


44 ±12 


2.2 ±0.6 


6.3 ±3.0 


4.8 ±2.2 




7.9 ±3.8 


M3 


15 ±8 


0.11 ±0.06 


2.2 ± 1.4 


1.6±1.1 




2.7±1.8 


M5 


15 ±7 


0.20 ±0.09 


2.2 ± 1.3 


1.6 ±1.0 




2.7±1.6 


M13 


16 ±7 


0.24 ±0.11 


2.3 ±1.3 


1.7 ±1.0 




2.9 ±1.7 


Ter 5 


100 ± 20 


2.6 ±0.5 


14 ±6 


11 ±5 




18 ±8 


NGC 6440 


68 ±30 


0.75 ± 0.33 


9.8 ±5.8 


7.3 ±4.3 




12 ±7 


NGC 6517 


30 ±15 


0.21 ±0.10 


4.3 ±2.7 


3.2 ±2.1 




5.4 ±3.4 


M28 


85 ±28 


2.0 ±0.7 


12 ±6 


9.2 ±4.7 




15.3 ±7.8 


NGC 6752 


27 ± 12 


1.1 ±0.5 


3.9 ±2.3 


2.9 ±1.7 




4.9 ±2.9 


M15 


53 ±20 


0.40 ±0.15 


7.6 ±4.1 


5.7±3.1 




9.5 ±5.2 




Model 3 


(minimum x^): 


^JL = -0.52 and 


<T = 0.68 Nio = 354 ± 43 




47 Tuc 


37 ±10 


2.5 ±0.7 


5.3 ±2.5 


4.0 ±1.9 




6.6 ±3.1 


MS 


12 ±6 


0.13 ±0.06 


1.7± 1.1 


1.3 ±0.8 




2.2 ±1.4 


M5 


13 ±6 


0.24 ±0.11 


1.9 ±1.1 


1.4 ±0.8 




2.3 ±1.4 


M13 


14 ±6 


0.30 ±0.13 


2.0 ±1.2 


1.5 ±0.9 




2.5 ±1.4 


Ter 5 


82 ±16 


2.9 ±0.6 


12 ±5 


8.9 ±3.8 




15 ±6 


NGC 6440 


48 ±21 


0.74 ± 0.33 


6.9 ±4.0 


5.2 ±3.0 




8.6 ±5.1 


NGC 6517 


23 ±12 


0.21 ±0.11 


3.3 ±2.1 


2.5 ±1.6 




4.1 ±2.7 


M28 


63 ±21 


1.5 ±0.5 


9.1 ±4.6 


6.8 ±3.5 




11 ±6 


NGC 6752 


21 ±10 


1.2 ±0.5 


3.0 ±1.8 


2.3 ±1.4 




3.8 ±2.3 


M15 


41 ±15 


0.43 ±0.16 


5.9 ±3.2 


4.4 ±2.4 




7.4 ±3.9 


.dio emission 




our 


study. For 47 Tuc, 


iMcConnell et al||2004|: 



A potentially very useful additional constraint comes from 
observations of the diffuse radio emission in GCs. Assuming 
that the only contribution to this flux is from the pulsars, 
then such measurements constrain the integrated luminosity 
function in a given cluster. Our Monte Carlo models make 
specific predictions for these observations (see Eqn. |4|. For 
Terzan 5, the total radio flux S'obs.tot = 5.2 mJy kpc^ (sum 
of diffuse flux a n d the fluxes of point sources) found by 
iFruchter fc GossI l|200d ). iFruchter fc Gosd l|2000l ) observed 
some other clusters too, among which NGC 6440 belongs 
to our list (Table |4)|. But as they mentioned that their ob- 
servation in this cluster is consistent with the position of a 
single pulsar PSR B1745— 20, we can not use this datum for 



Sobs, tot = 2.0 ± 0.3 mJy kpc . 

Assuming both the diffuse flux measurements for 
Terzan 5 and 47 Tuc are dominated by their respective pul- 
sar populations, we can confront them with the predictions 
from our simulations. An inspection of Tables [2] and [3] shows 
that the observed diffuse flux for 47 Tuc is successfully repro- 
duced by all of the models, within the nominal uncertainties. 
For Terzan 5, the power law models provide a better match 
to the diffuse flux overall, while the log-normal models pre- 



^ This sum is essentially equivalent to the individual fluxes of 
the 14 pulsars in this cluster with measured fluxes so far (see Ta- 
ble 1). The remaining 9 currently known pulsars, must therefore 
contribute much less than a mJy of diffuse flux. For example a 
typical flux of 30 ^Jy per pulsar would bring the diffuse flux to 
~ 2.3 mJy. 
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Table 3. Population estimations and predictions using power-law luminosity functions. Sec Table 2 for details about tabulated parameters. 



Cluster 



iV., 



(mjy) 



-'^7, Sim 

(r;^) = 0.08 (r;^> = 0.06 

(lO^"* erg s-i) 



Model 4: Lmin = 0.017 and ^ = 0.92 TVio = 3399 ± 421 



47 Tuc 
M3 
M5 
M13 
Ter 5 
NGC 6440 
NGC 6517 
M28 

NGC 6752 
M15 



313 ±84 
114 ±57 
112 ±50 
118 ±52 
764 ±153 
485 ± 217 
238 ± 119 
628 ± 209 
216 ±97 
411 ± 155 



3.4 ±0.89 
0.19 ±0.09 
0.32 ±0.14 
0.40 ±0.18 

4.4 ±0.9 
1.2 ±0.5 

0.35 ±0.18 

2.5 ±0.8 
1.9 ±0.9 

0.68 ±0.25 



45 ±21 
16 ±10 

16 ±10 

17 ±10 
110 ±48 
70 ±41 
34 ±22 
90 ±46 
31 ±19 
59 ±32 



34 ± 16 
12 ±8 

12 ±7 

13 ±8 
83 ±36 
53 ±31 
26 ±16 
68 ±35 
23 ±14 

449 ± 24 



Model 5: Lmin = 0.022 and /3 = 1.01 TVio = 3767 ± 478 



47 Tuc 
M3 
M5 
M13 
Ter 5 
NGC 6440 
NGC 6517 
M28 

NGC 6752 
M15 



324 ± 87 
121 ±61 
116 ±52 
123 ± 55 
815 ± 163 
580 ± 260 
271 ± 136 
714 ±238 
237 ± 106 
466 ±176 



3.3 ±0.9 

0.19 ±0.09 

0.31 ±0.14 

0.39 ±0.18 

4.5 ±0.9 

1.3 ±0.6 

0.38 ±0.20 

2.7 ±0.9 

2.0 ±0.9 

0.73 ±0.27 



47 ±22 
17 ±11 

17 ±10 

18 ±11 
117 ±51 
84 ±50 
39 ±25 
103 ± 53 
34 ±20 
67 ±36 



35 ± 17 
13 ±8 
13 ±7 
13 ±8 
88 ±39 
63 ±37 
29 ± 19 
77 ±30 
26 ± 15 
50 ±27 



Model 6 (FGOO): L^ 



0.03 and /3 = 0.85 A^io = 1590 ± 194 



47 Tuc 
M3 
M5 
M13 
Ter 5 

NGC 6440 
NGC 6517 
M28 

NGC 6752 
M15 



153 ± 41 
55 ±28 
54 ±24 
57 ±26 
363 ± 73 
216 ±96 
112 ±56 
287 ± 95 
101 ± 45 
192 ± 73 



3.2 ±0.86 

0.18 ±0.09 

0.31 ±0.14 

0.38 ±0.17 

4.2 ±0.8 

1.0 ±0.5 

0.33 ±0.16 

2.2 ±0.7 

1.8 ±0.8 

0.62 ±0.23 



22 ±10 
7.9 ±5.1 
7.8 ±4.6 
8.2 ±4.9 
52 ±23 
31 ±18 
16 ±10 
41 ±21 
15 ±9 
28 ±15 



17 ±8 
5.9 ±3.8 
5.8 ±3.4 
6.2±3.7 
39 ± 17 
23 ± 14 

12 ±8 
31 ± 16 

11 ±6 
21 ± 11 



(r,^> = 0.1 



56 ±27 
21 ±13 
21 ±12 
21 ±13 
138 ± 60 
87 ±52 
43 ±27 
113 ±58 
39 ±23 
74 ±40 



58 ±28 
22 ±14 

21 ±12 

22 ±13 
147 ±64 
104 ±62 
49 ±31 
129 ± 66 
43 ±25 
84 ±46 



28 ±13 
9.9 ±6.3 
9.7±5.7 
10 ±6 
65 ±29 
39 ±23 
21 ±13 
52 ±26 
18 ±11 
35 ±19 



100 



"*«SSj 



10 




10 



L [mJy kpc"] 



Figure 5. Observed luminosity distribution with simulated lu- 
minosity distributions generated with exponential distribution 
function with A = 0.676 mjy~^ kpc~-^ (lower curve) and A = 
0.439 mjy~^ kpc~^ (upper curve). "*"s represent the observed 
distribution. 



diet a slightly smaller flux that l i es 2-5 a below the nominal 
value found bv lFruchter fc Gosj i2000'). 

In Fig. [6] we fix /i to the nominal value from model 1, 
i.e. as found by FK06 for normal pulsars (-1.1) and vary 
a. For this case, we see that there are only two possible 
ranges of a which are compatible with the diffuse flux mea- 
surement of Terzan 5: o" ~ 0.5 or ct ~ 0.9 (see the upper 
panel of Fig. [6]). The "solution" with a ~ 0.5, however lies 
well outside the x^ contours shown in Fig. (2] We therefore 
favor the region with /x ~ —1.1 and a ~ 0.9 (the nominal 
FK06 values) which is consistent with both constraints. In 
this case, the implied total number of pulsars Ntot ~ 150 
(see the lower panel of Fig. [6|. Further constraints on these 
parameters using a more detailed Bayesian analysis of these 
constraints for Terzan 5 will be the subject of a subsequent 
paper (Chennamangalam et al. in preparation). 

Fig. [7] shows the analogous diagram to Fig. [6] for the 
power-law luminosity function for the choice /3 = 1. In this 
case, there is a wide range of Lmin values that are consistent 
with the diffuse flux measurements (upper panel of Fig. [7|), 
and no significant additional constraints on I/min can be 
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Table 4. Observed and inferred properties of globular clusters containing pulsars used in the present work. From left to right, we list the 
GC name, distance from the Sun d and from the Galactic center dgcon (both in kpc), concentration parameter c (base 10 logarithm of the 
ratio of the tidal radius to core radius), core radius r^ (pc), base-10 logarithm of the central density (pc in solar luminosities per cubic 
parsec), velocity dispersion Vc (km s~^), base-10 logarithms of the cluster mass {Mqq in solar masses) and base-10 logarithms of the 
core relaxation timescale (t(r'c) in yr), metallicity [F e/H], normalized two body encounter ra te Fnorm, gamma-ray flux L-y (10^^ erg s~^) 
and inferred number of gamma-ray pulsars A^.^ from lAbdo et all ll2010h and lTam et al.l 1 1201 j) . For Fnorm, we first calculate the two body 
encounter rate for each GC as F = p^'^ r^. The value s of core radi us have been calculated as re = d tan 9c where 9c are the angular radii 
as quoted in the latest version of the Harris catalog l|Harrijl 19961 . updated in December 2010). The central cluster density pc, has been 
calculated as pc = T,c/{rc p) using the values of central surface brightness Ijlvc (in V magnitude per square arcsecond) and the extinction 
coefficient Ay = 3.1 E(B — V) where E(B — V) is the color excess. Ec, the central surface brightness i n Lyp, pc^ can b e calculated 
as log(Sc) = 0.4 [26.392 - {pvc - Ay)] and p is a parameter defined as log(p) = -0.603 X 10"'= -I- 0.302 llDiorgovskilll993^ . Finally we 
normalized F to Fnorm considering F = 100 for M62 ffoUow inglAbdo et a [20ld). The tabulated values for Vc and log(MGc) can be found 
at http://www.astro.lsa. umich.edu/r^ognedin/gc/vesc.dat JGnedin et al.l 20021) . 



GC 



log(pc) Vc log(AfGc) log(t(rc)) [Fe/H] Fnor: 



N^ 



47Tuc 


4.03 


7.4 


2.07 


0.42 


4.93 


16.4 


6.17 


M3 


10.23 


12.0 


1.89 


1.10 


3.58 


9.2 


5.98 


M5 


7.76 


6.2 


1.73 


0.99 


3.87 


11.8 


5.93 


M13 


7.13 


8.4 


1.53 


1.28 


3.55 


10.3 


5.89 


Ter 5 


5.50 


1.2 


1.62 


0.25 


5.26 


12.7 


5.57 


NGC 6440 


8.47 


1.3 


1.62 


0.34 


5.23 


21.6 


5.91 


NGC 6517 


10.60 


4.2 


1.82 


0.19 


5.30 


20.6 


5.72 


M28 


5.70 


2.7 


1.67 


0.40 


4.86 


16.3 


5.74 


NGC 6752 


4.42 


5.2 


2.50 


0.22 


5.01 


7.1 


5.50 


M15 


10.30 


10.4 


2.29 


0.42 


5.08 


13.4 


6.08 



7.84 


-0.72 


41.27 


4.8ti:l 


33tl^ 


8.31 


-1.50 


2.61 


— 


— 


8.28 


-1.29 


5.83 


— 


— 


8.51 


-1.53 


3.26 


— 


— 


7.57 


-0.23 


46.50 


25.7tri 


180+J«« 


7.60 


-0.36 


78.27 


19.0+_f,,^ 


i3oir 


6.92 


-1.23 


29.71 


— 


— 


7.62 


-1.32 


28.19 


6.2t?:« 


43t^t 


6.88 


-1.54 


14.72 


i-4tS;? 


lot- 


7.84 


-2.37 


66.81 


<5.8 


<56 



made. The upper bound of S%ot gives an extremely high 
value of A'^tot ~ 1000 which seems unrealistic, but the lower 
bound of Stot gives A'tot = 340 for Lmin ~ 0.05 mjy kpc^ 
(lower panel of Fig. [7}. 



4.2 Predicted population sizes and diffuse fluxes 
for different GCs 

The detection of diffuse gamma-ray emission from GCs has 
allowed some constraints to be placed on A'^^, the number 
of g amma-ray emitting pulsars in each cluster. Following 
lAbd o ct aJi (,2010 ) . we can write the total gamma-ray lumi- 
nosity 



N,{E) iv^), 



(8) 



where (E) is the average spin-down power of MSPs, (ri-y) is 
the average spin-down to gamma-ray luminosity conversion 
efhciency. As the values of (E) and (n^) are not well known, 
I Abdo et all (|2010l ) assumed (E) = (1.8 ±0.7) x 10^* erg s'^ 
and (riy) = 0.08. We use the values of A'^^ estimated using 
the above relationship for the clusters with gamma-ray flux 
(and hence luminosity) measurements. 

In Fig. [SI we compare the estimates of A'^ with our pre- 
dicted numbers of radio pulsars (A^'rad) for one model of each 
luminosity function. A reasonable agreement can be noted 
for the log-normal function, but for the power-law function, 
the values of A^rad are significantly larger than those of A'^,. 
This fact remains unchanged even if we choose other mod- 
els from these luminosity functions (see Tables [5] and [3)). 
Although this simple analysis does provide some support 
to the log-normal models, due to the assumptions made in 
equation ((SJ and implicitly assuming that A'^, — A'lad, it 
does not help to constrain their values significantly. 



By assuming N^, = A^rad, we obtain L^, 



from equa- 



tion ([8]). We tabulate I/^.sim for different choices of {r]-y} for 



different models in Tables [5] and (3] These values can be com- 
pared with observed values of L^ and 7V^ as shown in Table 
31 It is apparent that the 7-ray luminosities predicted by 
the power-law models are generally higher than observed. 
We note, however, that in addition to the explicit assump- 
tions about beaming geometry mentioned above, it has been 
recently shown that the 7-ray observations can be biased by 
one or more very bright pulsars in the cluster (jFreire et al.l 
20111), and may not be representative of the diffuse flux of 
the whole population. 



5 COMPARISON WITH EARLIER RESULTS 

In Fig. [9] we compare our predicted number of pulsars hav- 
ing L > 0.5 mJy kpc^ in different GCs using FK06 param- 
eters to those by HCTIO. For Terzan 5, we do not use the 
distance they adopted for this cluster (10.3 kpc). Instead, 
here we re-calculate the value of N{L > 0.5) by adopting 
exactly the s ame method as HCT IO using the recent esti- 
mate 5.5 kpc (jOrtolani et aj|2007l ) to calculate luminosities. 
The overall agreement is good which again highlights the 
fact that, at least above 0.5 mJy kpc^, the exact form of the 
luminosity function for GC pulsars is not uniquely specified 
by the current sample of luminosities. 

HCTIO used their power-law luminosity functions to 
search for correlations between the number of inferred ra- 
dio pulsars and fundamental cluster parameters. In their 
Fig. 3, they present evidence for a correlation between A'^rad 
and both cluster metallicity [Fe/H] as well as the two-body 
encounter rate Fnorm. An inspection of these diagrams sug- 
gests that the claimed correlations are strongly influenced 
by Terzan 5. Adopting, for the purposes of this discussion, 
the parameters of model 1 (i.e. the log-normal luminosity 
parameters found by FK06), together with the revised dis- 
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L|,|jjj [mJy kpi:" 




L,,,,,, [lily tpt^^l 



Figure 6. Variations of total radio flux (upper panel) and pre- 
dicted number of pulsars (lower panel) in Terzan 5 as obtained 
from our simulations with log-normal distribution keeping fi fixed 
at —1.1 and varying a. The dashed lines in the upper panel denote 
the ranges of Ssi^.tot = •S'obs,tot ± 25%. 



Figure 7. Variations of total radio flux (upper panel) and pre- 
dicted number of pulsars (lower panel) in Terzan 5 as obtained 
from our simulations with Power law distribution keeping /3 fixed 
at 1.0 and varying Lmin- The dashed lines in the upper panel 
denote the ranges of Ssim.tot = 5obs,tot ± 25%. 



tance to Terzan 5, we revisit these proposed correlations in 
Fig. 1101 Also shown here are the results of correlation tests 
between A'rad and other cluster parameters. We searched 
for relationships between the distance of each GC and the 
galactic center, dgcon, the logarithm of the central luminosity 
density, pc, the concentration parametecl, c, the logarithm 
of the core relaxation time t{rc), the cluster mass, Mac, the 
central velocity dispersion, Vc, and the core radius, Vc- The 
parameter values used for this analysis are given in Table ID 

As can be seen in Fig. 1101 none of the scatter diagrams 
provides compelling evidence for a direct relationship be- 
tween A'rad and any cluster parameters. The lack of any sta- 
tistically significant correlations can also be seen formally 
in Table O where we have calculated Pearson's correlation 
coefficient (rp) and probability at which the null hypothesis 
of zero correlation is disproved (Prp), Spearman correlation 
coefficient (rs) and probability at which the null hypothesis 
of zero correlation is disproved (Prs), Kendall's r and the 
probability at which the null hypothesis of zero correlation 
is disproved {Pt}- 

We have shown the results of correlation analyses only 
for FK06 (our model 1), but for other models the results are 
almost the same. Because both Spearman correlation and 
Kendall's r test are based on ranks, and in all the models, 
the GCs with descending order of ranks (based on A^rad) 
are as follows Terzan 5, M28, NGC 6440, M15, 47 Tuc, 

^ This parameter is defined to be the logarithm of the ratio of 
the GCs tidal radius to its core radius. 



Table 5. Results of different statistical correlation tests between 
the predicted number of pulsars with various GC parameters. 



GC 


Pearson 


Spearman 


Kendall 


property 


rp 


Jr^rp 


rs 


Prs 


r 


Pr 


i norm 


0.60 


0.07 


0.80 


0.02 


0.67 


0.01 


Fe/H 


0.51 


0.13 


0.44 


0.18 


0.31 


0.21 


rfgccn 


-0.64 


0.04 


-0.69 


0.04 


-0.54 


0.03 


log{Pc) 


0.64 


0.04 


0.62 


0.06 


0.45 


0.07 


c 


-0.22 


0.55 


-0.24 


0.48 


-0.16 


0.52 


log(t(rc)) 


-0.29 


0.41 


-0.51 


0.13 


-0.34 


0.17 


Mac 


-0.25 


0.48 


-0.23 


0.49 


-0.22 


0.37 


Vc 


0.31 


0.38 


0.58 


0.08 


0.31 


0.21 


Tc 


-0.60 


0.07 


-0.52 


0.12 


-0.40 


0.10 



NGC 6517, NGC 6752, M13, M3 and M5 (with the ex- 
ception in model 1 and 2 where there is a tie between M3 
and M 5, but the order of other GCs are the same, see tables 
[3 and 131). Even for the parametric test - Pearson's correla- 
tion analysis, rp lies always within 12% of that of model 1 
for Fe/H and within 14% of that of model 1 for Fnorm- This 
also explains why our result contradicts with that of HCTIO 
inspite of overall good agreement between predicted num- 
ber - according to HCTIO, the GCs with descending order 
of ranks are as follows Terzan 5, 47 Tuc, M 28, NGC 6440, 
NGC 6441, NGC 6752, M13, M5, M3, which is diff'erent 
from what we obtain. 



© 0000 RAS, MNRAS 000, 000-000 



12 Bagchi, Larimer & Chennamangalam 



log(p^ in Lsun pc' ) 



c=log(r,/rJ 



, [kms-'l 



■"c [pi:] 



Mqc [ 1 0- Mjj,„] log(t(rJ in yr) ii|,„„ [kpc] 

Figure 10. Plot of predicted number of pulsars in different GCs using FK06 parameters against various cluster parameters (see text). 



6 CONCLUSIONS 

We have modeled the observed luminosity distribution of 
millisecond pulsars in globular clusters as the brighter tail 
of a parent distribution. We have found that either a log- 
normal or a power-law distribution can be used as the par- 
ent distribution. We have demonstrated that a wide range of 
possible luminosity functions are compatible with the data, 
and that log-normal distribution functions provide a bet- 
ter match to the data than the traditionally favored power- 
law distributions. In the light of these results, we conclude 
that there is currently no need to assume that the luminos- 
ity function for cluster pulsars is any different than that of 
pulsars in the Galactic disk found by FK06. Based on this 
result, it is quite possible that all pulsars follow similar lumi- 
nosity distribution irrespective of their positions or recycling 
history. 

Contrary to earlier claims by HCTIO, we find no evi- 
dence for a significant correlation between the inferred num- 
bers of radio pulsars in GCs and either metallicity or stel- 
lar encounter rate. No significant correlations were found 
against other cluster parameters either. Despite the lack of 
any obvious correlations found among this sample of 10 GCs, 
it is of great interest to perform an analysis using a much 
larger sample of clusters. Further constraints may be possi- 
ble by incorporating observations of the diffuse gamma-ray 



flux, though this approach is complicated by model depen- 
dencies in gamma-ray efficiency and the radio/gamma-ray 
beaming fraction. 

One key difference between the two luminosity functions 
we have not commented on thus far is shown in Tables[2]and 
[3] by the tabulated parameter A'^io, the sum of the popula- 
tion estimates across all 10 GCs. As can be seen, the power- 
law models predict a systematically larger parent population 
than for the log-normal distribution (i.e. A''io in the range 
1600-3800 compared to 350-700). This observation implies 
that the power-law distributions require larger birth rates 
over the log-normal models by a factor of 2-10. Although 
we shall defer a detailed population size analysis to a future 
paper, containing population estimates for more GCs, a sim- 
ple scaling of these numbers to all 150 GCs currently known 
implies a population range for potentially observable recy- 
cled pulsars of 5000-11000 pulsars in the log- normal models 
versus 24000-57000 for the power law models. Assuming a 
recycled pulsar lifetime of ~ 10^" yr, and a mean beaming 
fraction of 50%, the implied birth rate of this population is 
at least 10~^ yr~^ over all Galactic GCs. Recent results con- 
cerning the low-mass X-ray binary (LMXB) population in 
GCs (see lHeinkell201ll : |Poolevll2010l . for reviews) suggest that 
there are of order 200 LMXBs in Gala ctic GCs. Assuming a 
typical LMXB lifetime of order 10* yr (|Kulkarni k, Naravanl 
1 1988) ■ the implied birthrate is comparable to our rough es- 
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Figure 8. Comparison of predicted number of pulsars from our 
simulations witii those obtained from 7-ray fluxes. Top: log- 
normal luminosity function (FK06 parameters). Bottom: power- 
law luminosity function (Model 5). 
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Figure 9. Comparison of our predicted number of pulsars in 
different GCs using FK06 parameters to those by HCTIO. 



timates for the recycled pulsars, provided that the pulsar 
population estimates are closer to the ranges suggested by 
the log-normal models. 

We consider this study to be the first step towards a 
more comprehensive analysis of the pulsar content of GCs. 
More detailed studies of the pulsar luminosity functions 
which better account for the selection effects and detection 
issues in the various radio surveys are still needed to further 
probe all these issues. In particular, a search for correlations 
between cluster parameters and the pulsar content beyond 
the small sample of 10 GCs considered here is needed to 
better understand this diverse population of neutron stars. 
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